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We employed femtosecond time- and angle-resolved photoelectron spectroscopy to analyze the response 
of the electronic structure of the 122 Fe-pnictide parent compounds Ba / EuFe2As2 and optimally doped 
BaFei.85Coo.i5As2 near the T point to the optical excitation by an infrared femtosecond laser pulse. The 
experiments were carried out at an equilibrium temperature of 100 and 300 K with a pump fluence below 
1.4mJ/cm 2 . We identify pronounced changes of the electron population within several lOOmeV above 
and below the Fermi level, which we explain as combination of (i) coherent lattice vibrations, (ii) a hot 
electron and hole distribution, and (iii) transient modifications of the chemical potential. The response of the 
three different materials is very similar. In the Fourier transformation of the time-dependent photoemission 
intensity we identify three modes at 5.6, 3.3, and 2.6 THz. While the highest frequency mode is safely assigned 
to the Ai 9 mode, the other two modes require a discussion in comparison to literature. The time-dependent 
evolution of the hot electron distribution follows a simplified description of a transient three temperature 
model which considers two heat baths of lattice vibrations, which are more weakly and strongly coupled 
to transiently excited electron population. Still the energy transfer from electrons to the strongly coupled 
phonons results in a rather weak, momentum-averaged electron-phonon coupling quantified by values for 
A (a; 2 ) between 30 and 70meV 2 . The chemical potential is found to present pronounced transient changes 
reaching a maximum of 15meV about 0.6 ps after optical excitation and a damped oscillatory behavior 
resembling the intensity oscillations assigned to coherent phonons. This change in the chemical potential is 
particularly strong in a two band system like in the 122 Fc-pnictide compounds investigated here due to the 
pronounced variation of the electrons density of states close to the equilibrium chemical potential. 

PACS numbers: 74.70.Xa, 78.47.J-, 79.60.-i, 71.20.-b 



I. INTRODUCTION 



The discovery of high-T c superconductivity in F doped 
LaFcAsO [1\ led to an almost immediate discovery of fur- 
ther FeAs-based superconductors with transition temper- 
atures T c up to 55 K. While the superconducting pairing 
mechanism in conventional superconductors is related to 
electron-phonon coupling, current thinking is that the 
iron pnictides (FePn), similar to the cuprate supercon- 
ductors, belong to the class of unconventional supercon- 
ductors, in which the mechanism is possibly related to 
correlation effects and/or magnetic excitations @. Al- 
though there are various similarities between cupratcs 
and FcPns, e.g. the nearness in the phase diagram of an- 
tiferromagnetism and superconductivity, there are also 
differences. While in the cuprates the parent compounds 
are antiferromagnetic Mott-Hubbard insulators, in the 
FcPns the parent compounds are antiferromagnetic met- 
als which indicates that correlation effects in the latter 
systems are reduced. This view is supported by e.g. x-ray 
absorption spectroscopy on FcPns resulting in a modcr- 
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ate value of the on-site Coulomb repulsion of U « 1.5 e V 
1. 

On the other hand there is no real consensus on the 
mechanism of high-T c superconductivity in FcPns and 
on the role of electron-phonon coupling. Experiments on 
the isotope effect in these compounds exhibit conflict- 
ing results @. Theoretical calculations on the electron- 
phonon coupling constant A yielded for the non-magnetic 
and the antiferromagnetic system A « 0.2 and A < 0.35, 
respectively Optical pump-probe studies with fem- 

tosecond time-resolution analyzed A by determining the 
second moment of the Eliashbcrg electron-phonon cou- 
pling function A(w 2 ) through the rate of energy transfer 
from the optically excited electrons to phonons 0,0]. Al- 
though the reported results vary in detail due to different 
assumptions regarding ui the resulting electron-phonon 
coupling constant A is found to be 0.1 - 0.2. In our re- 
cent femtosecond time- and angle-resolved photoemission 
spectroscopy (trARPES) study in which ultrashort 
light pulses in the UV spectral range were used, we de- 
rived for the second moment of the Eliashberg electron- 
phonon coupling function, A(w 2 ) = 90±40meV 2 . Us- 
ing an average phonon frequency (uj) « 20meV/?i a 
value of A < 0.2 can be derived from this result. Al- 
though these electron-phonon coupling constants are 
small, which makes it difficult to explain T c values above 
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20 K by electron-phonon coupling, there is no consensus 
on the role of phonons in the mechanism of superconduc- 
tivity in FcPns. Therefore information on phonons and 
their coupling to the charge carriers is vitally important. 
In this view we present here further results on the coher- 
ent excitations of phonons and on the electron-phonon 
coupling in undoped Ba/EuFe2As2 and and Co doped 
Ba2As2 122 compounds using trARPES. 

Photoelectron spectroscopy is widely considered to be 
a surface sensitive method due to the small penetration 
depth of the free electron like final state into the material 
for energies of 10 - 100 eV. For the low photon energy 
of 6 - 7 eV used in laser photoemission experiments [8j- 
[Toj as well as in the present study a bulk sensitivity can 
be considered @ due to a larger mean free path of ex- 
cited electrons. However, on the one hand, a sensitivity 
to electronic surface states for these small kinetic ener- 
gies of photoelectrons is well established [ill; 13 ■ On the 
other hand, combined surface and bulk sensitivity 13 1 
as well as detection of states at buried interfaces 14 1 
were reported. Therefore, the degree of surface vs. bulk 
sensitivity might be highly material specific and the pho- 
toemission matrix elements need to be considered in de- 
tail The present study investigates 122 compounds 
with their pronounced two-dimensional electronic struc- 
ture. On the other hand 122 compounds exhibit an in- 
creasing three-dimesionality with increasing doping [l6j]. 
As will be presented below, the reported phenomena do 
not depend on doping. Therefore, the question of surface 
vs. bulk sensitivity appears not to be an essential one in 
the present context and for the following we consider to 
probe with trARPES the surface near region. 



II. EXPERIMENTAL DETAILS 

trARPES experiments were performed on EuFe2As2, 
BaFe2As2 parent compounds and an optimally doped 
BaFei.85Coo.i5As2 (T c = 23 K). Single crystals of 
EuFe2 As2 were grown by the Bridgman method [17J while 
Ba-based com pou nds were grown from self-flux in an alu- 
mina crucible |18| . All samples were cleaved in ultrahigh 
vacuum (p < 10~ 10 mbar) at T = 100 K where the major 
part of the measurements were carried out. Additional 
measurements on the parent compounds were performed 
at 300 K. 

A commercial regenerative amplifier Coherent RegA 
9050 generates ultrashort infrared laser pulses of 820 nm 
wavelength (1.5 eV) with a repetition rate of 300 kHz 
and pulse duration of 55 fs. The beam is then split in 
two parts, where one is used to provide the pump pulses 
and the other one is frequency-doubled and compressed 
twice to achieve 205 nm wavelength (6 eV) pulses with 
a duration of 80 fs, used for probing. We detect photo- 
electrons which are generated by the UV probe pulse 
using an electron time of flight spectrometer with an 
acceptance angle of ±3 degrees. Thereby we monitor 
the momentum and energy dependent single-particle 
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FIG. 1: (a) Color coded intensity plot of the trARPES intensity 
of BaFe1.85Coo.15 As2 at the T point as a function of binding energy 
and pump-probe delay, taken with an incident pump fluence of F = 
1.4 mj/cm 2 at T = 100 K. Red and green dashed lines mark the 
spectra shown in (b). (b) trARPES spectra for the first minimum 
(red) and maximum (green) of the oscillation in comparison to a 
spectrum before excitation (black). The solid lines are fits to the 
data (sec scction f!!! Bl) . The Fermi level position Ep extracted from 
the fit is indicated by the red and green arrows. 



spectral function in the vicinity of the Fermi level, 
Ep, as function of pump-probe delay. Measurements 
were carried out in normal emission, i.e. around the 
T-point. The energy resolution of 50meV is determined 
by the time of flight spectrometer and the bandwidth of 
the probe pulses. The overall temporal resolution was 
100 fs S. 



III. RESULTS AND DISCUSSION 

The time-resolved photoelectron intensity of 
BaFci.85Coo.i5As2 at the r-point is shown in figure 
[IJa) as function of binding energy and time delay for 
an incident excitation fluence of F = 1.4mJ/cm 2 . 
After optical excitation, (i) pronounced oscillations 
modulate the photoemission spectrum around Ep and 
(ii) optically excited electrons and holes broaden the 
distribution function around Ep. In this report we 
disentangle these contributions and discuss the coherent 
phonon and electron dynamics in sections IIII Al and 
IIII Bl respectively. Figure db) compares spectra taken 
at At — 100 fs and At = 200 fs, corresponding to the 
first minimum and maximum of the oscillations, with 
a spectrum taken before excitation. Apart from the 
strong change of the distribution function at Ep after 
excitation, a rigid shift of the cutoff at Ep induced by 
the oscillation is observed. 



A. Coherent phonon dynamics 

We start with the discussion of the pronounced oscil- 
lations recognized in figure QJa) . The time-dependent 
photoemission intensity integrated for energies E > Ep, 
I{t), is depicted in figure (Ha) for different pumping flu- 
ences. This integration yields a rather good signal-to 
noise ratio due to the high contrast at the Fermi cutoff. 
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FIG. 2: (a) Integrated trARPES intensity of BaFei.85Coo.isAs2 
as a function of pump-probe delay for four different fiuenccs mea- 
sured at 100 K. Dashed lines are the fitted spline function used 
to obtain background subtraction. In panel (b) the respective 
background-subtracted oscillations are shown. 



Within the first picosecond after optical excitation, an 
oscillation with a period close to 200 fs is observed. For 
At > 1.5 ps, an interference pattern is found which indi- 
cates the presence of further modes. This is particularly 
evident for the data at higher fluence. These oscillations 
are superimposed on a steep rise in integrated intensity 
after the optical excitation followed by an intensity de- 
crease. In the following we refer to the oscillatory part as 
the coherent contribution and to the non-oscillatory part 
as the incoherent contribution to I(t). To determine the 
coherent contribution, the incoherent contribution was 
subtracted from I(t). As a first step, the incoherent 
part was described by a single exponential decay. While 
this description was adequate at low excitation densi- 
ties, for higher excitation densities this procedure does 
not provide a satisfying background determination. For 
this reason, we adopted a method consisting of fitting a 
smooth spline function to the data, which was originally 
developed for the analysis of extended x-ray absorption 
fine structure (EXAFS) [l9| and was recently used to ex- 
tract the coherent oscillations of the charge density wave 
(CDW) amplitude mode in TbTe 3 H3- This method 
provides a smooth background while leaving the oscil- 
lation unaffected as shown in figure EJb) . For more de- 
tails see [2(| ■ After background substraction, fast Fourier 
transformation (FFT) was used to obtain the frequency 
spectrum of the coherent contribution AI co h(t). Figure 
[3Ja) depicts the FFTs for the corresponding fluences pre- 
senting a sharp intense peak at cj-y = 5.6 THz (23±1 meV) 
and two weaker modes at CJ2 = 3.3 THz (14 ± 1 meV) and 
ujs = 2.6 THz (11 ± meV). As discussed below these fea- 
tures are likely to be fingerprints of coherent lattice vibra- 
tions and are observed in trARPES due to the coupling 
between the electronic system and the coherent phonons. 
All three modes were found with similar frequencies also 
in the parent compounds BaFe2As2 and EuFe 2 As2 (fig- 
ure [S^b)) suggesting that these coherent excitations are 
a general phenomenon in 122 FePns. 

The resulting frequencies are summarized in figure|4]as 
a function of the incident pumping fluence for the three 



samples under investigation. Open and closed symbols 
denote data taken at T — 300 K and T = 100 K, respec- 
tively. The experimental errors are determined by the 
point spacing of the FFT which depends on the extend 
of the respective data set in the time domain. Additional 
uncertainties arise from the decay of the oscillations in 
the time domain, that leads to the broadening of peaks 
in the FFT. We find that the determined FFT frequen- 
cies do not vary with the excitation density in all three 
cases and are temperature-independent (see also figure 
|3jc)). The FFT amplitude for uj\ in BaFei.85Coo.i5As 2 
in the lower panel of figure 2] is well described by a lin- 
ear dependence with fluence. In combination with the 
fluence independent frequencies, this points to the linear 
regime of the coherent excitations, showing no indication 
of anharmonicity within our sensitivity 21 1. 

To determine the geometry and character of the coher- 
ent modes, we compare our data with Raman scattering 
results [22rl27l | and a theoretical study Q. By this com- 
parison, we clearly identify the mode uj\ at 23 meV with 
the Raman active As Ai g phonon mode, corresponding 
to a displacement of the As atoms perpendicular to the 
FeAs layers, as sketched in the inset of figure [3ta). This 
mode was found in Raman scattering studies of various 
pnictide 122 compounds [221 - I271 ] . including BaFe2As2. In 
addition, coherent excitation of this mode has also been 
observed in time-resolved reflectivity measurements on 
Co doped BaFe2As2 [|| and in terahertz spectroscopy on 
BaFe2As2 (28|. The frequency was found to be 5.56 THz 
in the first case and 5.5 THz in the latter, which is in 
agreement with our observation. A more recent study 
found a coherent excitation of the Ai ff mode also in Sm- 
FeAsO [29|], which suggests that the coherent excitation 
of this mode is a general phenomenon in the FePns. How- 
ever, the other two modes were not observed with time- 
resolved optical spectroscopy. 

For the excitation of coherent phonons, traditionally 
two models have been established [30l - l32l . the displacive 
excitation of coherent phonons (DECP) [33[ and the im- 
pulsive stimulated Raman scattering (ISRS) [34|. While 
the DECP mechanism is supposed to be the main exci- 
tation mechanism in opaque media like metals and semi- 
metals and allows to first order only the excitation of 
fully symmetric Ai modes [33| . the ISRS is commonly 
accepted as the excitation process in transparent media 
like insulators and some semiconductors and allows for 
the excitation of all Raman active modes @, Hl| . It was 
shown that ISRS can also be extended to opaque ma- 
terials and is the more general description [30(. From 
the observation of coherently excited phonons here in 
trARPES, we can conclude on a strong coupling of the 
coherent vibration to electronic states directly at Ep, as 
we measure the imprint of the coherent modes on the 
electronic system. 

The assignment of the other two modes is less clear. A 
possible candidate for the UJ2 mode could be the E 9 mode, 
a shear vibration of the Fe and As atoms within the FeAs 
plane (26j . On the one hand this Raman active mode was 
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FIG. 3: (a) Fast Fourier amplitudes of the data in figure [2] for correspondent flucnecs showing the three coherent modes oji=5.6 THz, 
UJ2=3. 3 THz, u}3=2.6 THz. The inset is a sketch of the Aig mode, corresponding to a vertical displacement of the As atoms with respect 
to the Fe plane, (b) Comparison of FFTs for BaFei.g5Coo.i5As2, BaFe2As2 and EuFe2As2 at T = fOO K. Traces are vertically offset for 
clarity, (c) Comparison of FFTs of BaFe2As2 for T=300 K and T=f00 K. All three modes are found to be independent of T. 
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FIG. 4: Upper panel: FFT frequencies of the three modes 
as function of excitation fluence F for EuFe2As2 (triangles), 
BaFei.85Coo.i5As2 (circles) and BaFe2As2 (diamonds). Open and 
closed symbols represent data taken at T = 300 K and T = f 00 K, 
respectively. Error bars are defined by the inter-point spacing of 
the FFT. Horizontal lines are guide to the eyes. Within error bars, 
the frequencies of the modes do not depend on the excitation den- 
sity. Lower panel: FFT amplitude for U)i of BaFei.85Coo.i5As2 as 
function of incident pumping fluence. The dashed line is a linear 
fit to the data. 



found at 16meV |26( corresponding to 3.8 THz, which is 
considerably higher than the mode ui2 in our experiment 
and this frequency for the E 9 mode is consistent with Ra- 
man observations of this mode in BaFe2 As2 [HI ■ On the 
other hand, a Raman scattering study on CaFe2As2 (HJ 
reported the observation of a mode at 3.35 THz, which 
could be attributed to the E ff mode. In addition, a recent 
study of Co doped BaFe2As2 using inelastic x-ray scat- 
tering found modes close to the BZ center with an energy 
s=s 13meV [3a|. These values are all consistent with our 



observation. 

Comparing our frequencies with the calculated phonon 
dispersion in we see that another possible candidate 
for the LU2 mode could be the mode at the X-point around 
13 meV, which shows combined in-plane and out-of-plane 
oscillations with wave vector q equal to the nesting wave 
vector Q n and enhanced electron-phonon (e-ph) coupling 
according to the calculations. The coherent excitation of 
phonon modes with finite momenta is usually prohibited, 
as the exciting photons allow only vertical transitions due 
to their negligible momentum. However, a specific mo- 
mentum dependent optical excitation probability could 
in principle generate a transient momentum dependent 
population of hot charge carriers that could drive also 
bosons modes with finite momenta. Furthermore, con- 
tributions of surface phonons or spin excitations [33, [3c| , 
which do not show up in the calculations or in Raman 
spectroscopy, cannot be ruled out. 

Regarding the w 3 mode, the calculation in Q shows 
a mode at the zone-center around 10 meV which is, 
however, not Raman active. Furthermore, according 
to the calculation, this mode couples only weakly to 
the electronic system and thus would not be observed 
here. Therefore additional work is required in order to 
unambiguously clarify the geometry and character of 
these two modes, u>2 and uj^. 

We will now concentrate on the time evolution of the 
three modes. Considering the FFTs in figure^ we notice 
that the mode uj\ exhibits a considerably broader peak in 
the frequency domain than the other two modes UJ2 and 
u>3. This indicates a faster relaxation in the time domain 
of this mode. In order to study this time evolution in 
more detail, we concentrated on the BaFei.85Con.i5As2 
data using the highest pump fluence of F = 1.4mJ/cm 2 
and applied a sliding window fast Fourier transform 
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FIG. 5: (a) Color coded intensity plot of the swFFT intensity of 
BaFei.85Coo.i5As2 (at F = 1.4 mj/cm 2 , T = 100 K) as function 
of frequency and time delay of the sliding window center. A sliding 
window width of 2.4 ps was used, (b) swFFT corresponding to four 
different sliding window centers, marked by arrows in (a), (c) Time 
evolution of the integrated swFFT amplitude of u\ and LJ2+LO3 
integrated together, plotted on a logarithmic scale. Solid lines are 
fits to the data. 



(swFFT) analysis, with a time-window of 2.4 ps which 
allows a frequency resolution of ~ 0.5 THz, sufficient to 
resolve the two modes L02 and W3 that are closest to each 
other. The intensity of the swFFT signal as function of 
the transformed frequency and the sliding window center 
is shown in figure [5ta). No change in the frequency of 
the modes with time is observed, which is consistent with 
the harmonic regime of the oscillations concluded above 
and from figure |U In addition, the amplitude of u>\ is 
indeed damped with a higher decay rate than the am- 
plitude of the other two modes. This becomes obvious 
in figure E^b), where the swFFT intensity correspond- 
ing to four different sliding window centers, marked by 
arrows in figure [5Ja) , is shown. The integrated swFFT 
amplitude at the frequencies of the modes is plotted on 
a logarithmic scale as a function of the window position 
in figure 02c). Note that due to interference effects be- 
tween u>2 and UJ3 , combined amplitudes corresponding to 
these two modes were integrated. The integrated swFFT 
amplitude of both lo\ and L02 + UJ3 shows a clear linear 
behavior in the logarithmic plot. By fitting single expo- 
nential decay functions to the respective amplitude time 
evolution (solid lines), we determine t± = 1.2 ± 0.1 ps 
and Tu, 2 +u 3 = 3.1 ± 0.4 ps as decay times. We will now 
discuss the different relaxation times for the mode ui\ 
and the modes W2 and 0J3. The relaxation of coher- 
ently excited phonons occurs both through pure dephas- 
ing and population decay of the excited modes. While 
the pure dephasing takes place mainly by scattering at 
crystal defects and e-ph interaction, the population de- 



cay is governed by anharmonic phonon-phonon (ph-ph) 
interaction [3§] . The scattering with crystal defects is un- 
likely to explain the different decay behavior of different 
modes within the same crystal, assuming a homogeneous 
excitation density of these modes. This could be differ- 
ent for e.g. surface modes. The other two relaxation 
processes can be distinguished by their temperature and 
fluence dependence. On the one hand, the relaxation 
by anharmonic ph-ph decay into two acoustic phonons 
shows an increase in decay rate with increasing temper- 
ature [ll|, [3!|. This is not observed in our case, as we 
have not found an increase of the width of the peaks in 
the FFT spectra for the higher temperature as shown 
in figure G3c). On the other hand, the scattering of co- 
herent phonons with excited carriers by e-ph coupling is 
highly dependent on the density of excited carriers and 
thus on the excitation density, which leads to a faster 
decay at higher excitation densities [2l|, |39| . This is in- 
deed observed in the data in figurc[3]Ja), where we find an 
increase with fluence of the linewidth of oj\ . Therefore, 
we conclude that the dominating relaxation channel for 
the coherent oscillations is the coupling to excited carri- 
ers by e-ph scattering. The shorter lifetime of the mode 
loi compared to the other two modes thus indicates a 
stronger coupling to the electronic system for the mode 
lo\. This is consistent with the lower amplitude of the 
other two modes, as e-ph coupling is the key parame- 
ter for both the generation and the detection of coherent 
phonons in trARPES and, as a consequence, the ampli- 
tude of coherent phonons in the detected signal scales 
with the coupling strength. 

Interestingly, the relaxation that we have found for the 
mode lo\ is clearly faster than in the time-resolved opti- 
cal reflectivity measurements, where the Ai s mode was 
found to relax with a time-constant TA lg ~ 2.5 ps 
This might be explained by a difference in crystal quality, 
exhibiting less crystal defects in the optical pump-probe 
experiments. Another possibility might be the different 
detection technique used. Whereas the optical reflectiv- 
ity probes the entire region of the optical penetration 
depth of several tens of nanometers, photoelectron spec- 
troscopy is a surface sensitive technique. Thus, the co- 
herent phonons in the sub-surface region might exhibit a 
faster relaxation than in the bulk of the lattice. 



B. Dynamics of the electronic system 

So far we have analyzed the oscillations in time- 
dependent photoelectron intensity. As it has been men- 
tioned above, two further phenomena are observed in the 
data presented in figurc[T] These are (i) the photo-excited 
distribution of electrons and holes and (ii) a variation of 
the high-energy cutoff of the spectra. While the distri- 
butions of hot charge carriers have been discussed be- 
fore already in detail for other materials j40l - l44l | the clear 
change in the high energy cutoff is less established. Here 
we report the results of an analysis based on the assump- 
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FIG. 6: Timc-dcpcndent electronic temperature T e (a) and chem- 
ical potential fi (b) determined by fitting the transient trARPES 
spectra shown in figure^ (see text). The results of a fit with a 
three temperature model (3TM) for T e (red), the temperature of 
a heat bath of more strongly coupled phonon modes T p (dashed 
green) and the lattice temperature 2j (dashed-dotted dark violet) 
are also shown in (a). 



tion that these two phenomena can be described by the 
combination of a thermalized hot electron distribution 
function and a transient modification of the chemical po- 
tential. The fits presented in figure fflh) are the result 
of this model approach and describe the data very well. 
These fits have been performed as a function of time de- 
lay and the resulting hot electron temperature T e (i) and 
chemical potential fx(t) are depicted in figure |Hta) and 
EKb), respectively. Note that for this analysis a constant 
auxiliary density of states was assumed and non-thermal 
hot electrons which occur before the thermalization of the 
electronic system close to to were neglected here for sim- 
plicity [55|. We find a strong modulation of fi(t) which 
indicates that the intensity oscillations discussed above 
originate from a time-dependent change fi = fi(t) and 
that the coherent phonon excitation modifies the chem- 
ical potential. Although details of this phenomenon re- 
main to be analyzed in future work we tentatively ex- 
plain the changes in /z(t) as a consequence of the modu- 
lation of the electronic structure in response to coherent 
phonons [111, [HI, [TH, [4f| [46| . Due to the significant varia- 
tion of the electronic density of states near Ep in the two 
band structure of FePns systems [13], charge neutrality 
might well require a transient readjustment of fj, to the 
instantaneous position of the electron and hole pocket at 
X and r, respectively. 

To describe the time-dependent change of T e we em- 
ploy a three temperature model (3TM) as used be- 
fore for the High T c superconductor Bi2Sr2CaCu20s+5 
(BSCCO) (4l[. This description gives a reasonable agree- 
ment with the obtained T e (t) as shown in figure [SJa). 



A similar analysis within the same 3TM has been re- 
cently performed for FePns using timc-resolvcd optical 
reflectivity measurements [48[. The 3TM describes the 
preferential coupling of the optically excited hot electron 
distribution with temperature T e to a more strongly cou- 
pled subset of phonon modes with temperature T p , which 
subsequently couple on a longer timescale to the rest of 
the lattice modes with temperature 7] . The following set 
of coupled differential equations is used [49| : 



at 

dTp 

at 

m 

dt 



-H(T e ,T p ) + — 

-9±H(T e ,T p )-^i 
Cp Tp — Ti 



(1) 

(2) 
(3) 



The source term S is a Gaussian excitation pulse with 
a FWHM of 55 fs and an energy density determined by 
the absorbed fluence distributed over the optical penetra- 
tion depth. The specific heat of electrons is determined 
by C e = 7 e T e , with the electronic specific heat coefficient 
7 e . For the specific heat of the hot phonons and the 
rest of the phonons, the Einstein model with a mode ujq 
is used, C p = Sfhuo fr| Tp and Q = 3(1 - f)hu f^, 

where n = (e Bu, °/ fcBT — is the Bose-Einstein distribu- 
tion function. The parameter / is the fraction of modes 
that are more strongly coupled to the electrons. The 
anharmonic decay (responsible for the loss of oscillation 
signal) of the respective subset of hot phonons to the re- 
maining, more weakly coupled phonons, is described by 
Tp, For the energy transfer from the electrons to the 
more strongly coupled phonons, the formula derived by 
Allen [H is used: 



H{T e , Tp) 



3h\ (lj 2 ) T e - T p 



nkf 



T, 



(4) 



where A (w 2 ) is the second moment of the Eliashbcrg 
e-ph coupling function a 2 F(io). The c-ph coupling to the 
weakly coupled fraction (1 — /) barely influences the evo- 
lution of T e and has therefore been neglected. Electronic 
diffusion processes can be assumed not to be essential 
here as the electric and thermal conductivity is strongly 
reduced along the longer c-axis due to the layered, quasi 
2D structure of the FePns [10, HH - 

Naturally, the description of the phonon system with 
a simple Einstein mode is a strong simplification and 
a proper choice of the mode energy influences the re- 
sult of the simulation. To provide an estimate for the 
frequency uiq, we consider the phonon spectrum of the 
122 FeAs compounds. The calculated phonon spectrum 
extends up to ~ 35meV as derived in several calcula- 
tions [EIISII3 an d supported by experiments [47L l52l l53l| 
with the strongest contributions to a 2 F(uj) between 
I0-30meV [|. 
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The 3TM was fitted to the transient T e for various val- 
ues of luo between 10 — 30meV. The resulting values for 
the parameters of the model are summarized in Table U 
for the three investigated compounds. For / and A (w 2 ), 
the left and right values are for a choice of ojq = lOmeV 
and 30 meV, respectively. An exemplary fit of the elec- 
tronic temperature T e using luq = 18meV is shown in 
figure [5] (a) for the Co doped sample as solid line, which 
shows excellent agreement with the measured data. For 
other choices of ujq , similar good agreement to the data is 
found. The temperatures of hot phonons (T p ) and other 
phonons (T;) derived from the model are shown as dashed 
and dash-dotted lines, respectively. 



compound 


/ A 


y) 


(meV 2 ) 


Tp (ps) 


EuFe2As2 


0.43 - 0.38 


56 


-65 


2.9 ±0.3 


BaFei.85Coo.isAs 2 


0.48 - 0.43 


40 


-55 


3.3 ±0.5 


BaFe 2 As 2 


~ 0.45 


30 


-46 


~ 3.5 



TABLE I: Fitting parameters of the 3TM. Left and right values for 
/ and A correspond to a choice of ajo = 10 meV and 30 meV, 

respectively. 



The values determined for A (w 2 ) can be compared to 
recently published values for the e-ph coupling strength 
obtained from optical pump-probe experiments. Mansart 
ct al. [48| report values of A^cj 2 ) w 64meV 2 for almost 
optimally doped BaFe2-zCo x As2,which was determined 
by the 3TM. Stojchevska et al. [jj derived a somewhat 



higher value of A (lj 2 ) = 110 ± lOmeV 2 for SrFe2As2 
from the temperature dependence of quasiparticle relax- 
ation times and for SmFeAsOi-^F^ an even larger value 
of A(w 2 ) = 135 ± lOmeV 2 is reported [Hj]. In general, 
our results are compatible with the weak e-ph coupling 
found in the FePns and indicate even weaker coupling for 
BaFe2As2 compounds than in EUFC2AS2. However, one 
has to bear in mind that the performed analysis assumes 
a thermalized electronic system and systematically ne- 
glects non-thermal electrons. As thermalization of non- 
thermal hot electrons leads to an additional heating of 
the thermalized electron distribution, this process can 
counteract the energy transfer to the lattice and leads to 
a slower energy relaxation of the thermalized electrons. 
Thus, the analysis of only the thermalized part of the 
electronic distribution might considerably underestimate 
e-ph coupling especially at early times. 

The fraction / « 0.4 of preferentially coupled modes is 
in agreement with the observations of Mansart et al. [48| . 
and is considerably higher than found in BSCCO by 
/ 



0.2 41 



Here we follow the explanation given by 
Mansart et al., who suggest that this stronger selectivity 
of e-ph coupling in the cuprates can be explained by the 
stronger 2D character of the cuprates, which could favor 
the selective coupling to specific modes. The FePns and 
especially Co doped 122 compounds are characterized 
by a stronger interlayer coupling and a more 3D char- 



pling 



acter 54 1, which might lead to a less selective e-ph cou- 



48[ . The equilibration time Tp of the hot phonon 



distribution with the rest of the lattice modes is consider- 
ably faster here than in the optical reflectivity measure- 
ments by Mansart et al., who report Tp ~ 5 — 7ps. This 
faster relaxation of excited phonon modes observed in 
trARPES measurements compared to optical reflectivity 
is consistent with the faster relaxation of the coherently 
excited phonon modes discussed earlier and indicates a 
different phonon relaxation behavior in the bulk and the 
surface/subsurface region, respectively. In addition, the 
excitation densities in the experiments by Mansart et al. 
where a factor of two to three larger than here, and an 
increase of Tp with fluence was observed (48j . 

Based on our results of A {lj 2 ) we can estimate the 
value of the e-ph coupling constant A for a particular 
value of lj. Considering the coherently excited Ai g mode 
at 23 meV, which showed enhanced coupling in calcula- 
tions, we find A < 0.15 for all compounds. This estimate 
is in agreement with calculations jl, Q of various FePns 
compounds, which report average values of A < 0.35. 
Taking the mean of the phonon DOS as reference, A gets 
even smaller in agreement with other publications. Even 
if we consider the lowest coupled modes around 12meV 
to be most important for e-ph coupling, A does not ex- 
ceed a value of 0.5. Similar small values for A have been 
found in the cuprates [4lj | , which suggests limited impor- 
tance of e-ph coupling for the pairing mechanism in both 
classes of materials. 



IV. CONCLUSION 

Femtosecond time- and angle-resolved photoelectron 
spectroscopy was demonstrated to be a powerful tool to 
investigate the response of the electronic system in the 
vicinity of the Fermi level. In this study we could sepa- 
rate for various 122 FePns compounds transient changes 
in the electron distribution function and the chemical 
potential. The first one was analyzed by a simplified 
three temperature model. In agreement with earlier 
time-resolved optical studies we conclude a rather small 
momentum-averaged electron phonon coupling constant 
A < 0.15, which suggests a limited importance of e-ph 
coupling for the Cooper pair formation in these systems. 
The transient chemical potential is explained to be the 
consequence of charge neutrality in the probed region 
of the material while the electronic system is in a non- 
equilibrium state characterized by hot electrons and co- 
herent as well as incoherent phonons. Such changes in the 
chemical potential can be expected to be a general phe- 
nomena, however, we consider that the two band nature 
of the 122 FePns compounds enhances the effect in com- 
parison to a single band system. The pronounced effect 
of coherent phonons on the electronic structure near Ep 
demonstrate vivid electron-phonon coupling but the as- 
signment of the two observed lower frequency excitations 
beside the Ai s mode requires further attention. Future 
studies as function of doping might provide more insight 
into the origin of the changes in the chemical potential 



8 



as well as the material specificity of all observed coherent knowledges the Alexander von Humboldt Foundation, 
phonon modes. We thank Hermann Diirr and Lilia Bocri for fruitful dis- 

cussions. 

ACKNOWLEDGEMENTS 



Funding from the Deutsche Forschungsgemeinschaft 
within SPP 1458 is gratefully acknowledged. R.C Ac- 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono. 
Iron-based layered superconductor LaOi-^F^FeAs (a; = 
0.05-0.12) withT c = 26K. J. Am. Chem. Soc, 130:3296, 
2008. 

[2] G. R. Stewart. Superconductivity in iron compounds. 
Rev. Mod. Phys., 83:1589-1652, 2011. 

[3] T. Kroll, S. Bonhommeau, T. Kachel, H. A. Duerr, 
J. Werner, G Behr, A. Koitzsch, R. Huebel, S. Leger, 
R. Schoenfelder, A. K. Ariffin, R. Manzke, F. M. F. 
de Groot, J. Fink, H. Eschrig, B. Buechner, and 
M. Knupfer. Electronic structure of LaFeAsOi-zF^ from 
x-ray absorption spectroscopy. Phys. Rev. B, 78:220502, 
2008. 

[4] L. Boeri, O. V. Dolgov, and A. A. Golubov. Is 
LaFeAsOi-ajFa: an electron-phonon superconductor? 
Phys. Rev. Lett, 101:026403, 2008. 

[5] L. Boeri, M. Calandra, I. I. Mazin, O. V. Dolgov, and 
F. Mauri. Effects of magnetism and doping on the 
electron-phonon coupling in BaFe2As2. Phys. Rev. B, 
82:020506, 2010. 

[6] B. Mansart, D. Boschetto, A. Savoia, F. Rullier- 
Albenque, A. Forget, D. Colson, A. Rousse, and 
M. Marsi. Observation of a coherent optical phonon in 
the iron pnictide superconductor Ba(Fe 1 _ a .Coa;) 2 As2 {x 
= 0.06 and 0.08). Phys. Rev. B, 80:172504, 2009. 

[7] L. Stojchevska, P. Kusar, T. Mertelj, V. V. Kabanov, 
X. Lin, G H. Cao, Z. A. Xu, and D. Mihailovic. Electron- 
phonon coupling and the charge gap of spin-density wave 
iron-pnictide materials from quasiparticle relaxation dy- 
namics. Phys. Rev. B, 82:012505, 2010. 

[8] L. Rettig, R. Cortes, S. Thirupathaiah, P. Gegenwart, 
H. S. Jeevan, M. Wolf, J. Fink, and U. Bovensiepen. 
Ultrafast momentum-dependent response of electrons in 
antiferromagnetic EuFe2As2 driven by optical excitation. 
Phys. Rev. Lett, 108:097002, 2012. 

[9] T. Kiss, F. Kanetaka, T. Yokoya, T. Shimojima, 
K. Kanai, S. Shin, Y. Onuki, T. Togashi, C. Zhang, C. T. 
Chen, and S. Watanabe. Photoemission spectroscopic ev- 
idence of gap anisotropy in an /-electron superconductor. 
Phys. Rev. Lett, 94:057001, 2005. 
[10] J. D. Koralek, J. F. Douglas, N. C. Plumb, Z. Sun, 
A. Fedorov, M. M. Murnane, H. C. Kapteyn, S. T. Cun- 
diff, Y. Aiura, K. Oka, H. Eisaki, and D. S. Dessau. 
Laser based angle-resolved photoemission, the sudden 
approximation, and quasiparticle-like spectral peaks in 
BiaSraCaCuaOs+A". Phys. Rev. Lett, 96:017005, 2006. 
[11] P. A. Loukakos, M. Lisowski, G. Bihlmayer, S. Bliigel, 
M. Wolf, and U. Bovensiepen. Dynamics of the self- 
energy of the Gd(0001) surface state probed by fem- 
tosecond photoemission spectroscopy. Phys. Rev. Lett., 



98:097401, 2007. 

[12] Uwe Bovensiepen, Hrvoje Petek, and Martin Wolf. Dy- 
namics in Solid States Surfaces and Interfaces, Vol. 1, 
Current Developments. Wiley- VCH, Weinheim, 2010. 

[13] L. Perfetti, P. A. Loukakos, M. Lisowski, U. Bovensiepen, 
H. Berger, S. Biermann, P. S. Cornaglia, A. Georges, and 
M. Wolf. Time evolution of the electronic structure of 
lT-TaS2 through the insulator-metal transition. Phys. 
Rev. Lett, 97:067402, 2006. 

[14] L Rettig, P S Kirchmann, and U Bovensiepen. Ultrafast 
dynamics of occupied quantum well states in pb/si(lll). 
New J. Phys., 14:023047, 2012. 

[15] W. Schattke. Prog. Surf. Scl, 64:89, 2000. 

[16] S. Thirupathaiah, S. de Jong, R. Ovsyannikov, H. A. 
Diirr, A. Varykhalov, R. Follath, Y. Huang, R. Huisman, 
M. S. Golden, Yu-Zhong Zhang, H. O. Jeschke, R. Va- 
lentf, A. Erb, A. Gloskovskii, and J. Fink. Orbital charac- 
ter variation of the Fermi surface and doping dependent 
changes of the dimensionality in BaFea-iCo^Asa from 
angle-resolved photoemission spectroscopy. Phys. Rev. 

B, 81:104512, 2010. 

[17] H. S. Jeevan, Z. Hossain, Deepa Kasinathan, H. Rosner, 

C. Geibel, and P. Gegenwart. Electrical resistivity and 
specific heat of single-crystalline EuFe2As2: A magnetic 
homologue of SrFe 2 As 2 . Phys. Rev. B, 78:052502, Aug 
2008. 

[18] F. Hardy, P. Adelmann, T. Wolf, H. v. Lohneysen, and 
C. Meingast. Large anisotropic uniaxial pressure depen- 
dencies of T c in single crystalline Ba(Fe g2 Coo.o8) 2 As2. 
Phys. Rev. Lett, 102:187004, 2009. 

[19] P. A. Lee, P. H. Citrin, P. Eisenberger, and B. M. 
Kincaid. Extended x-ray absorption fine structure~its 
strengths and limitations as a structural tool. Rev. Mod. 
Phys., 53:769-806, Oct 1981. 

[20] F. Schmitt, P.S. Kirchmann, U. Bovensiepen, R. G. 
Moore, J-H Chu, D. H. Lu, L. Rettig, M. Wolf, I. R. 
Fisher, and Shen Z-X. Ultrafast electron dynamics in 
the charge density wave material TbTe3. New J. Phys., 
13(063022):26, 2011. 

[21] M. Hase, M. Kitajima, S. Nakashima, and K. Mizoguchi. 
Dynamics of coherent anharmonic phonons in bismuth 
using high density photoexcitation. Phys. Rev. Lett., 
88:067401, 2002. 

[22] L. Chauviere, Y. Gallais, M. Cazayous, A. Sacuto, M. A. 
Measson, D. Colson, and A. Forget. Doping dependence 
of the lattice dynamics in Ba(Fe 1 _ a ,Co a: ) 2 As2 studied by 
Raman spectroscopy. Phys. Rev. B, 80:094504, 2009. 

[23] L. Chauviere, Y. Gallais, M. Cazayous, M. A. Measson, 
A. Sacuto, D. Colson, and A. Forget. Raman scatter- 
ing study of spin-density- wave order and electron-phonon 



9 



coupling in Ba(Fe 1 _ x Co x ) 2 As 2 . Phys. Rev. B, 84:104508, 
2011. 

[24] K.-Y. Choi, D. Wulferding, P. Lemmens, N. Ni, S. L. 
Bud'ko, and P. C. Canfield. Lattice and electronic 
anomalies of CaFe2As2 studied by Raman spectroscopy. 
Phys. Rev. B, 78:212503, 2008. 

[25] P. Kumar, A. Bera, D. V. S. Muthu, A. Kumar, U. V. 
Waghmare, L. Harnagea, C. Hess, S. Wurmehl, S. Singh, 
B. Buechner, and A. K. Sood. Raman evidence for the 
superconducting gap and spinphonon coupling in the su- 
perconductor Ca(Fe Q 95 Coo.o5) 2 As2. J. Phys.: Condens. 
Matter, 23:255403, 2011. 

[26] A. P. Litvinchuk, V. G. Hadjiev, M. N. Iliev, Bing Lv, 
A. M. Guloy, and C. W. Chu. Raman-scattering study of 
KsSri-zFeaAsa (x = 0.0, 0.4). Phys. Rev. B, 78:060503, 
2008. 

[27] S. Sugai, Y. Mizuno, K. Kiho, M. Nakajima, C. H. Lee, 
A. Iyo, H. Eisaki, and S. Uchida. Erratum: Pairing 
symmetry of the multiorbital pnictide superconductor 
BaFei.84Coo.i6As2 from raman scattering. Phys. Rev. B, 
83:019903, 2011. 

[28] K. W. Kim, A. Pashkin, H. Schaefer, M. Beyer, 
M. Porer, T. Wolf, C. Bernhard, J. Demsar, R. Hu- 
ber, and A. Leitenstorfer. Ultrafast transient genera- 
tion of spin-density-wave order in the normal state of 
BaFe2As2 driven by coherent lattice vibrations. Nat. 
Mat, doi:10.1038/nmat3294, 2012. 

[29] T. Mertelj, P. Kusar, V. V. Kabanov, L. Stojchevska, 
N. D. Zhigadlo, S. Katrych, Z. Bukowski, J. Karpin- 
ski, S. Weyeneth, and D. Mihailovic. Quasiparticle re- 
laxation dynamics in spin-density-wave and supercon- 
ducting SmFeAsOi-zFa; single crystals. Phys. Rev. B, 
81:224504, 2010. 

[30] G. A. Garrett, T. F. Albrecht, J. F. Whitakcr, and 
R. Merlin. Coherent thz phonons driven by light pulses 
and the sb problem: What is the mechanism? Phys. Rev. 
Lett, 77:3661, 1996. 

[31] M. Hase, K. Ishioka, J. Demsar, K. Ushida, and M. Ki- 
tajima. Ultrafast dynamics of coherent optical phonons 
and nonequilibrium electrons in transition metals. Phys. 
Rev. B, 71:184301, 2005. 

[32] K. Ishioka & O. V. Misochko. Coherent lattice oscilla- 
tions in solids and their optical control. Progress in Ultra- 
fast Intense Laser Science, in Springer Series in Chemi- 
cal Physics., 98, 2010. 

[33] H. J. Zeiger, J. Vidal, T. K. Cheng, E. P. Ippen, G. Dres- 
selhaus, and M. S. Dresselhaus. Theory for displacive ex- 
citation of coherent phonons. Phys. Rev. B, 45:768-778, 
1992. 

[34] Y. X. Yan, E. B. Gamble, and K. A. Nelson. Impulsive 
stimulated scattering: General importance in femtosec- 
ond laser pulse interactions with matter, and spectro- 
scopic applications. J. Chem. Phys., 83:5391, 1985. 

[35] A. P. Litvinchuk, Bing Lv, and C. W. Chu. Raman scat- 
tering study of electron-doped Pr x Cai_ :E Fe2As2 super- 
conductors. Phys. Rev. B, 84:092504, 2011. 

[36] D. Reznik, K. Lokshin, D. C. Mitchell, D. Parshall, 
W. Dmowski, D. Lamago, R. Heid, K.-P. Bohnen, A. S. 
Sefat, M. A. McGuire, B. C. Sales, D. G. Mandrus, 
A. Subedi, D. J. Singh, A. Alatas, M. H. Upton, A. H. 
Said, A. Cunsolo, Yu. Shvyd'ko, and T. Egami. Phonons 
in doped and undoped BaFe2As2 investigated by inelastic 
x-ray scattering. Phys. Rev. B, 80:214534, 2009. 

[37] A. Melnikov, A. Povolotskiy, and U. Bovensiepen. 



Magnon-enhanced phonon damping at Gd(0001) and 
Tb(0001) surfaces using femtosecond time-resolved op- 
tical second-harmonic generation. Phys. Rev. Lett., 
100:247401, 2008. 
[38] D. S. Inosov, J. T. Park, P. Bourges, D.L. Sun, Y. Sidis, 

A. Schneidewind, K. Hradil, D. Haug, C. T. Lin, 

B. Keimer, and V. Hinkov. Normal-state spin dynam- 
ics and temperature-dependent spin-resonance energy in 
optimally doped BaFei.g5Coo.i5As2- Nat Phys., 6:178- 
181, 2010. 

[39] M. Hase and Kitajima. M. Interaction of coherent 
phonons with defects and elementary excitations. Jour- 
nal of Physics: Condensed Matter, 22:073201, 2012. 

[40] M. Lisowski, P.A. Loukakos, U. Bovensiepen, J. Stahler, 

C. Gahl, and M. Wolf. Ultra-fast dynamics of electron 
thermalization, cooling and transport effects in Ru(001). 
Appl. Phys. A, 78(2), 2004. 

[41] L. Perfetti, P. A. Loukakos, M. Lisowski, U. Bovensiepen, 
H. Eisaki, and M. Wolf. Ultrafast electron relaxation 
in superconducting Bi2Sr2CaCu2 08+5 by time-resolved 
photoelectron spectroscopy. Phys. Rev. Lett., 99:197001, 
2007. 

[42] W. S. Fann, R. Storz, H. W. K. Tom, and J. Bokor. 
Electron thermalization in gold. Phys. Rev. B, 46:13592, 
1992. 

[43] N. Del Fatti, C. Voisin, M. Achermann, S. Tzortzakis, 

D. Christofilos, and F. Vallee. Nonequilibrium electron 
dynamics in noble metals. Phys. Rev. B, 61:16956-16966, 
2000. 

[44] H.-S. Rhie, H. A. Diirr, and W. Eberhardt. Femtosecond 
electron and spin dynamics in Ni/W(110) films. Phys. 
Rev. Lett, 90:247201, 2003. 

[45] F. Schmitt, P. S. Kirchmann, U. Bovensiepen, R. G. 
Moore, L. Rettig, M. Krenz, J.-H. Chu, N. Ru, L. Per- 
fetti, D. H. Lu, M. Wolf, I. R. Fisher, and Z.-X. Shen. 
Effect of the amplitude mode and the transient melting 
of a charge density wave on the electronic structure of 
tbte 3 . Science, 321(5896) :1649, 2008. 

[46] F. Schmitt, P. S. Kirchmann, U. Bovensiepen, R. G. 
Moore, J-H Chu, D. H. Lu, L. Rettig, M. Wolf, I. R. 
Fisher, and Z-X Shen. Ultrafast electron dynamics in 
the charge density wave material tbte3. New J. Phys., 
13:063022, 2011. 

[47] M. Zbiri, H. Schober, M. R. Johnson, S. Rols, R. Mit- 
tal, Y. Su, M. Rotter, and D. Johrendt. Ab initio lat- 
tice dynamics simulations and inelastic neutron scatter- 
ing spectra for studying phonons in BaFe2As2: Effect 
of structural phase transition, structural relaxation, and 
magnetic ordering. Phys. Rev. B, 79:064511, 2009. 

[48] B. Mansart, D. Boschetto, A. Savoia, F. Rullier- 
Albenque, F. Bouquet, E. Papalazarou, A. Forget, 
D. Colson, A. Rousse, and M. Marsi. Ultrafast tran- 
sient response and electron-phonon coupling in the iron- 
pnictide superconductor Ba(Fe 1 _ a ,Co :E ) 2 As2. Phys. Rev. 
B, 82:024513, 2010. 

[49] P. B. Allen. Theory of thermal relaxation of electrons in 
metals. Phys. Rev. Lett, 59:1460-1463, 1987. 

[50] M. A. Tanatar, N. Ni, C. Martin, R. T. Gordon, H. Kim, 
V. G. Kogan, G. D. Samolyuk, S. L. Bud'ko, P. C. Can- 
field, and R. Prozorov. Anisotropy of the iron pnictide su- 
perconductor Ba(Fe 1 _ ;c Co ;c ) 2 As2 (x = 0.074, T c = 23K). 
Phys. Rev. B, 79:094507, 2009. 

[51] J.-Ph. Reid, M. A. Tanatar, X. G. Luo, H. Shak- 
eripour, N. Doiron-Leyraud, N. Ni, S. L. Bud'ko, P. C. 



10 



Canfield, R. Prozorov, and Louis Taillefer. Nodes in 
the gap structure of the iron arsenide superconductor 
Ba(Fe 1 _ ;c Co a ;) 2 As2 from c-axis heat transport measure- 
ments. Phys. Rev. B, 82:064501, 2010. 

[52] S. E. Hahn, Y. Lee, N. Ni, P. C. Canfield, A. I. Gold- 
man, R. J. McQueeney, B. N. Harmon, A. Alatas, B. M. 
Leu, E. E. Alp, D. Y. Chung, I. S. Todorov, and 
M. G. Kanatzidis. Influence of magnetism on phonons 
in CaFe2As2 as seen via inelastic x-ray scattering. Phys. 
Rev. B, 79:220511, 2009. 

[53] R. Mittal, L. Pintschovius, D. Lamago, R. Heid, K-P. 
Bohnen, D. Reznik, S. L. Chaplot, Y. Su, N. Kumar, 
S. K. Dhar, A. Thamizhavel, and Th. Brueckel. Mea- 



surement of Anomalous Phonon Dispersion of CaFe2As2 
Single Crystals Using Inelastic Neutron Scattering. Phys. 
Rev. Lett, 102:217001, 2009. 

[54] P. Vilmercati, A. Fedorov, I. Vobornik, U. Manju, 
G. Panaccione, A. Goldoni, A. S. Sefat, M. A. McGuire, 
B. C. Sales, R. Jin, D. Mandrus, D. J. Singh, and N. Man- 
nella. Evidence for three-dimensional fermi-surface topol- 
ogy of the layered electron-doped iron superconductor 
Ba(Fe 1 _ ;r Co :c ) 2 As2. Phys. Rev. B, 79:220503, 2009. 

[55] For several systems we have observed and partly included 
such a subsystem of non-thermalized electrons [40l . l4lj . 



